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Abstract. Rapid uniform rotation of newborn neutron 
stars (protoneutron stars) is studied for a range of inter- 
nal temperatures and entropies per baryon predicted by 
the existing numerical simulations. Calculations are per- 
formed using general relativistic equations of hydrostatic 
equilibrium of rotating, axially symmetric stars. Stability 
of rotating configurations with respect to mass shedding 
and the axially symmetric perturbations is studied. Nu- 
merical calculations are performed for a realistic dense 
matter equation of state, under various assumptions con- 
cerning neutron star interior (large trapped lepton num- 
ber, no trapped lepton number, isentropic, isothermal). 
For configurations with baryon mass well below the max- 
imum one for the non-rotating models, the mass shed- 
ding limit depends quite sensitively on the position of 
the "neutrinosphere" (which has a deformed, spheroidal 
shape); this dependence weakens with increasing baryon 
mass. The absolute upper limit on rotation frequency is, to 
a good approximation, obtained for the maximum baryon 
mass of rotating configurations. Empirical formula for the 
maximum rotation frequency of uniformly rotating pro- 
toneutron stars is shown to be quite precise; it actually co- 
incides with that used for cold neutron stars. Evolutionary 
sequences at fixed baryon mass and angular momentum, 
which correspond to evolution of protoneutron stars into 
cold neutron stars are studied, and resulting constraints 
on the maximum rotation frequency of solitary pulsars are 
discussed. 
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1. Introduction 

Neutrons stars are born in gravitational collapse of mas- 
sive, degenerate stellar cores. Newly born neutron stars 
are hot and lepton rich objects, quite different from ordi- 
nary low temperature, lepton poor neutron stars. In view 
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of these differences, newly born neutron stars are called 
protoneutron stars; they transform into standard neutron 
stars on a timescale of the order of ten seconds, needed for 
the loss of a significant lepton number excess via emission 
of neutrinos trapped in the dense, hot interior. 

In view of the fact that the typical evolution timescale 
of a protoneutron star (seconds) is some three orders 
of magnitude longer, than the dynamical timescale for 
this objects (milliseconds), one can study its evolution 
in the quasistatic approximation (Burrows & Lattimcr 
1986). Properties of static (non-rotating) protoneutron 
stars, under various assumptions concerning composition 
and equation of state (EOS) of hot, dense stellar interior 
were studied by numerous authors (Burrows & Lattimcr 
1986, Takatsuka 1995, Bombaci et al. 1995, Bombaci 1996, 
Bombaci et al. 1996). 

The scenario of transformation of a protoneutron star 
into a neutron star could be strongly influenced by a phase 
transition in the central region of the star. Brown and 
Bethe (1994) suggested a phase transition implied by the 

condensation at supranuclear densities. Such a 
condensation could dramatically soften the equation of 
state of dense matter, leading to a low maximum allowable 
mass of neutron stars. In such a case, the massive pro- 
toneutron stars could be stabilized by the effects of high 
temperature and of the presence of trapped neutrinos, and 
this would lead to maximum baryon mass of protoneu- 
tron star larger by some 0.2 Mq than that of cold neu- 
tron stars. The deleptonization and cooling of protoneu- 
tron stars of baryon mass exceeding the maximum allow- 
able baryon mass for neutron stars, would then inevitably 
lead to their collapse into black holes. The dynamics of 
such a process was recently studied by Baumgarte et al. 
(1996). It should be mentioned, however, that the very 
possibility of existence of the kaon condensate (or other 
exotic phases of matter, such as the pion condensate, or 
the quark matter) at neutron star densities is far from 
being established. Recently, for instance, Pandharipande 
et al. (1995) pointed out that kaon-nucleon and nucleon- 
nucleon correlations in dense matter raise significantly the 
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threshold density for kaon condensation, possibly to the 

densities higher than those characteristic of stable neu- 
tron stars. In view of these uncertainties, we will restrict 
in the present paper to a standard model of dense matter, 
composed of nucleons and Icptons. 

The calculations of the static models of protoneutron 
stars should be considered as a first step in the studies of 
these objects. It is clear, in view of the dynamical scenario 
of their formation, that protoneutron stars are far from be- 
ing static. Due to the nonzero initial angular momentum 
of the collapsing core, protoneutron stars arc expected 
to rotate. On the other hand, the formation scenario in- 
volves compression (with overshoot of central density) and 
a hydrodynamical bounce, so that a newborn protoneu- 
tron star begins its life in a highly excited state, pulsating 
around its quasistatic equilibrium. In the present paper 
we study the rotation of protoneutron stars; pulsations of 
protoneutron stars will be discussed in a separate paper 
(Gondek, Haensel & Zdunik, in preparation). 

Some aspects of rapid uniform rotation of protoneu- 
tron stars have been recently studied in (Takatsuka 1995, 
Hashimoto et al. 1995). However, the calculations reported 
by Takatsuka (1995) were actually done for static (non- 
rotating) protoneutron stars, and were then used to esti- 
mate the maximum rotation frequency of uniformly rotat- 
ing protoneutron stars, fimax, via an "empirical formula". 
It should be stressed, that the validity of such an "em- 
pirical formula", which expresses fimax in terms of the 
mass and radius of the extremal static configuration with 
maximum allowable mass, had been checked only in the 
restricted case of cold neutron stars (Haensel &: Zdunik 
1989, Friedman et al. 1989, Shapiro et al. 1989, Haensel 
et al. 1995, Nozawa et al. 1996). Only isentropic equations 
of state were considered by Takatsuka (1995). Hashimoto 
et al. (1995) calculated the structure of stationary config- 
urations of uniformly rotating protoneutron stars, using a 
two-dimensional general relativistic code. These authors 
restricted themselves to the case with zero trapped lepton 
number. They assumed a constant temperature in the hot 
interior of the star, and used a zero temperature (cold) 
EOS for p < 10^" g cm^'^. It should be stressed, that the 
assumption of T = const, corresponds to an isothermal 
state in the Newtonian (flat space-time) theory of gravi- 
tation. In general relativity, we will define isothermal state 
by T* = = const, (where N is the lapse function and 
r is the Lorentz factor, see Section 3.1), and the effects of 
the space-time curvature will turn out to be rather impor- 
tant for massive neutron stars. Also, their choice for the 
low density edge of the hot interior can be questioned. Fi- 
nally, their criterion for finding maximally rotating config- 
uration is actually valid only for cold (T = 0) or isentropic 
protoneutron stars: its use in the case of the T = const. 
hot interior is unjustified (see Section 3.2 for a correct 
statement of the stability criterion). In a recent paper, 
Lai and Shapiro (1995) have studied the secular evolu- 
tion, secular "bar instability" , and the gravitational wave 



emission from the newly formed, rapidly rotating neutron 

stars. However, these authors used unrealistic (polytropic) 
equations of state of neutron star matter. Moreover, the 
calculations were done within Newtonian theory of gravi- 
tation. In view of this, the internal structure of their mod- 
els of newly born neutron stars was quite different from 
that characteristic of the realistic models of protoneu- 
tron stars. The problem of the secular "bar instability" 
in rapidly rotating neutron stars was also studied, using 
general relativity, by Bonazzola et al (1995). However, nu- 
merical calculations were done only for realistic equations 
of state of cold neutron star matter. 

In the present paper we study the properties of uni- 
formly rotating protoneutron stars, using exact relativis- 
tic description of the rapid, stationary rotation, combined 
with realistic equations of state of hot dense matter, used 
in the whole range of temperatures and densities rele- 
vant for protoneutron stars. In particular, we calculate 
the maximum frequency of uniform rotation of protoneu- 
tron stars and its dependence on their baryon mass, and 
on the thermal state and composition of stellar interior. 
It is clear, that uniform rotation represents only an ap- 
proximation to the actual rotational state of a newly born 
protoneutron star. Existing numerical simulations of grav- 
itational collapse of rotating cores of massive stars pro- 
duce differentially rotating protoneutron stars (Janka & 
Moenchmeyer 1989a, b, Moenchmeyer & Mueller 1989). 
However, it should be stressed that the initial rotational 
state of collapsing core is unknown, and this implies un- 
certainty concerning the rotational state of resulting pro- 
toneutron star. It is reasonable to say, that the actual 
degree of nonuniformity of rotation of a protoneutron star 
should be considered as unknown. In the present paper 
we will not address the question of the physical mecha- 
nisms that could "rigidify" the rotational motion within 
the protoneutron star interior. However, we will use the 
approximation of uniform rotation in order to limit the 
number the parameter space for our numerical calculation, 
and also because of the relative simplicity of the stability 
analysis in this specific, idealized case. 

Within our simplified model, the "ncutrinosphcrc" 
(which has actually a deformed, spheroidal shape) will 
separate hot, neutrino-opaque interior of a protoneutron 
star (hereafter referred to as "hot interior" ) from a signif- 
icantly cooler, neutrino-transparent envelope. The actual 
thermal state of the hot interior of protoneutron star is 
determined by its formation scenario, and is expected to 
be influenced by the dissipative processes (damping of pul- 
sations, viscous damping of differential rotation, neutrino 
diffusion). For simplicity, we will restrict ourselves to two 
limiting cases: an isothermal (T* = = const., see Sec- 
tion 3), and an isentropic (entropy per baryon s = const.) 
hot interior. We will also consider two limiting cases of 
the lepton composition of the protoneutron star interior. 
The flrst case, referring to the very initial state of pro- 
toneutron star, will correspond to a fixed trapped lepton 
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number. In the second case, neutrinos will not contribute 

to the Icpton number of the matter, which will correspond 
to vanishing chemical potential of the electron neutrinos; 
such a situation will take place after a deleptonization of 
a protononeutron star. The position of the neutrinosphere 
will be located using a simple prescription based on spe- 
cific properties of the neutrino opacity of hot dense matter. 
In all cases, the equation of state of hot dense matter will 
be determined using one of the models of Lattimer and 
Swesty (1991). 

The plan of the paper is as follows. In Section 2 wc 
describe the physical state of the interior of protoneutron 
star, with particular emphasis on the EOS of the hot in- 
terior at various stages of evolution. We explain also our 
prescription for locating the "neutrinosphere" of a pro- 
toneutron star, and we give some details concerning the 
assumed temperature profile within a protoneutron star. 
Using simple estimates of the timescales relevant for var- 
ious transport processes, we justify the approximation of 
stationarity which is used throughout this paper. In Sec- 
tion 3 we give a brief description of the exact equations, 
used for the calculation of stationary configuration of uni- 
formly rotating protoneutron stars. We discuss also stabil- 
ity of rotating configurations with respect to the axially- 
symmetric perturbations. The numerical method, used for 
the calculation of rapidly rotating configurations of pro- 
toneutron star, is briefly described in Section 4, where we 
also discuss numerical precision of our solutions. Maxi- 
mum rotation frequency, for various physical conditions 
prevailing in the hot stellar interior, calculated as a func- 
tion of the baryon (rest) mass of protoneutron star, is 
presented in Section 5. Then, in Section 6 we show the va- 
lidity of an empirical formula, which enables one to express 
with a surprisingly high precision the maximum frequency 
of rotating protoneutron stars in terms of the mass and 
radius of the maximum mass configuration of static (non- 
rotating) protoneutron stars with same EOS. In Section 
7 we study the evolutionary transformation of a rotat- 
ing protoneutron star into a cold neutron star. We show 
that, at fixed rest mass and angular momentum, maximum 
rotation frequency of protoneutron stars imposes severe 
constraints on the rotation frequency of solitary neutron 
stars. Finally, Section 8 contains discussion of our results 
and conclusion. 

2. Physical state of the interior of protoneutron 
stars 

We consider a protoneutron star (PNS) just after its 

formation. We assume it has a well defined "neutri- 
nosphere" , which separates a hot, neutrino-opaque interior 
from colder, neutrino-transparent outer envelope. Impor- 
tant parameters, which determine the local state of the 
matter in the hot interior are: baryon (nucleon) number 
density n, net electron fraction Yg = {rig- — nf,+ )/n, and 
the net electron-neutrino fraction = Y^^ — Y^^ . The cal- 
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culation of the composition of hot matter and of its EOS 
is described below. 

2.1. Neutrino opaque core with trapped lepton number 

Such a situation is characteristic of the very initial stage of 
existence of a PNS. Matter is composed of nucleons (both 
free and bound in nuclei) and leptons (electrons and neu- 
trinos; for simplicity, we do not include muons). All con- 
stituents of the matter (plus photons) are in thermody- 
namic equilibrium at given values of n, T and Yi = Y^+Yi,. 
The composition of the matter is calculated from the con- 
dition of beta equilibrium, combined with the condition of 
a fixed Yi, 

Yi = Ye + Y,, (1) 

where /zj are the chemical potentials of matter con- 
stituents. At the very initial stage we expect V; ~ 0.4. 
Electron neutrinos are degenerate, with fi^^ » T (in what 
follows we measure T in energy units). The deleptoniza- 
tion, implying the decrease of Yi , occurs due to diffusion 
of neutrinos outward (driven by the ^i,^ gradient), on a 
timescalc of seconds (Sawyer & Soni 1979, Bombaci ct al. 
1996). The diffusion of highly degenerate neutrinos from 
the central core is a dissipative process, resulting in a sig- 
nificant heating of the neutrino-opaque core (Burrows & 
Lattimer 1986). 

2.2. Neutrino opaque core with Y^ =0 

This is the limiting case, reached after complete delep- 
tonization. There is no trapped Icpton number, so that 
Yi =Ye and Y^^ = Yp^, and therefore iJ.,y^ = jio^ = 0. Neu- 
trinos trapped within the hot interior do not influence the 
beta equilibrium of nucleons, electrons and positrons, and 
for given n and T the equilibrium value of Yg is determined 
from 

IJ,p+ He = , (2) 

while fie+ = —fie- In practice, this approximation can be 
used as soon as electron neutrinos become non-degenerate 
within the opaque core, /i,/^, < T, which occurs after 
some ~ 10 seconds (Sawyer & Soni 1979, Bombaci et 
al. 1996). The neutrino diffusion is then driven by the 
temperature gradient, and the corresponding timescale of 
the heat transport (PNS cooling) can be estimated as 
~ 50(T/10 MeV)-'' s (Sawyer & Soni 1979). 

2.3. Neutrinosphere and temperature profile 

In principle, the temperature (or entropy per nucleon) 
profile within a PNS has to be determined via evolution- 
ary calculation, starting from some initial state, and tak- 
ing into account relevant transport processes in the PNS 
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interior, as well as neutrino emission from PNS. Trans- 
port processes within neutrino-opaque interior occur on 
timescales of seconds, some three orders of magnitude 
longer than dynamical timescales. The very outer layer 
of PNS becomes rapidly transparent to neutrinos, delep- 
tonizes, and cools on a very short timescale via e~e^ pair 
annihilation and plasmon decay. It seemed thus natural 
to model the thermal structure of the PNS interior by a 
hot core limited by a "neutrinosphere" , and an outer enve- 
lope of r < 0.5 MeV. The transition through the "neutri- 
nosphere" is accompanied by a temperature drop, which 
takes place over some interval of density just above the 
"edge" of the hot neutrino-opaque core, situated at some 

In view of the uncertainties in the actual tempera- 
ture profiles within the hot interior of PNS, we consid- 
ered two extremal situations for n > n^, corresponding 
to an isentropic and an isothermal hot interior. In the 
first case, hot interior was characterized by a constant 
entropy per baryon s = const.. In the case of trapped 
lepton number, this leads to the EOS of the type: pres- 
sure p = p{n, [s,Yi]), energy density e = e(n, [s,!!]), 
and temperature T = T{n, [s,Yi]), with fixed s and Y/. 
This EOS will be denoted by EOS[s, F/]. In the case of an 
isentropic, zero trapped lepton number EOS, we will have 

EOS[s,r^ = 0]. 

The condition of isothermality, which corresponds to 
thermal equilibrium, is more complicated. Due to the cur- 
vature of the space-time within PNS, the condition of 
isothermality corresponds to the constancy of T* = -^T 
(see Section 3.1). The significance of the T* ~ const. 
condition will be discussed in Section 3.1. In the static 
case, the isothermal state within the hot interior will be 
reached on a timescale corresponding to thermal equilibra- 
tion, which is much longer than the lifetime of a PNS. In 
the case of a rotating PNS, the situation can be expected 
to be even more complicated (see, e.g., Chapter 8 of Tas- 
soul (1978)). Nevertheless, we considered the T* = const. 
models for several reasons. First, as a limiting case so dif- 
ferent from the s ~ const, one, it enables us to check 
the dependence of our results for rapidly rotating PNS 
on the thermal state of the hot interior. Moreover, for 
the isothermal PNS we can apply the criterion of stability 
with respect to the axi-symmetric perturbations, and this 
will enable us to calculate the value of ilmax for the sta- 
ble supramassive rotating PNS models with an isothermal 
interior (see Section 3). 

Our calculation of the "neutrinosphere" within the hot 
PNS interior was done using a simple method, described 
below. For a given static PNS model, the neutrinosphere 
radius, R,y, has been located through the condition 

/ X / F ^ Vlg^dr = 1 , (3) 

where \u is the neutrino mean free path which is calcu- 
lated at the matter temperature T, while E^, is the mean 



energy of non-degenerate neutrinos at (and above) the 
neutrinosphere, — 3.15T^. We assumed that opacity 
above Ri, is dominated by the elastic scattering off nuclei 
and nucleons, so that = A°(n,T)/iJ^. Then, we de- 
termined the value of the density at the neutrinosphere, 
n^, for a given static PNS model, combining Eq. (||) with 
that of hydrostatic equilibrium, and readjusting accord- 
ingly the temperature profile within the outer layers of 
PNS. This value of n^ was then used in the calculations of 
rotating PNS models. Let us notice, that similar approxi- 
mation, in which the "neutrinosphere" in rapidly rotating 
protoneutron star was defined as a surface of constant 
density, was used by Janka and Moenchmeyer (1989a, b). 
Some details concerning actual calculation of the temper- 
ature profile in the vicinity of the "neutrinosphere" will 
be given in subsection 3.3. 

2.^. Equation of state and static models of PNS 

The starting point for the construction of our EOS for 
the PNS models was the model of hot dense matter of 
Lattimer and Swesty (1991), hereafter referred to as LS. 
Actually, we used one specific LS model, corresponding to 
the incompressibility modulus at the saturation density of 
symmetric nuclear matter K = 220 MeV. For n > n^ we 
supplemented the LS model with contributions resulting 
from the presence of trapped neutrinos of three fiavours 
(electronic, muonic and tauonic) and of the corresponding 
antineutrinos. 

In Fig. la, b we show our EOS in several cases, 
corresponding to various physical conditions in the hot, 
neutrino-opaque interior of PNS. For the sake of compar- 
ison, we have shown also the EOS for cold catalyzed mat- 
ter, used for the calculation of the (cold) NS models. In 
Fig. la we show EOS at subnuclear densities. At these den- 
sities, both the temperature and the presence of trapped 
neutrinos stiffen the EOS, as compared to the cold cat- 
alyzed matter one, and the stiffening is rather dramatic. 
The constant T EOS stiffens considerably at lower densi- 
ties, which is due to the weak dependence of the thermal 
contribution (photons, neutrinos) on the baryon density of 
the matter (this effect will be to some extent moderated 
by the factor in the isothermal PNS, see Section 3.1). It 
is quite obvious, that T = const. EOS becomes dominated 
by thermal effects below for n < 10^^ fm^^. On the con- 
trary, for isentropic EOS, the effect of the trapped lepton 
number (Yi = 0.4) turns out to be much more important 
than the thermal effects. This can be see in Fig. la, by 
comparing dash-dotted curve, [s — 2, Yi = 0.4], with the 
dotted line, which corresponds to an artificial (unphysical) 
case with small thermal effects, [s — 0.5, Y/ — 0.4]. 

It is clear, that the correct location of the "neutri- 
nosphere" , which separates hot interior from the colder 
outer envelope, should be important for the determination 
of the radius of PNS, and in consequence, of the maximum 
rotation frequency of a given PNS model. 
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Fig. 1. a. Pressure versus baryon density for our model of 
dense hot matter, under various physical conditions, for the 
subnuclear densities (n < 0.16 fm^^). The curve T = cor- 
responds to cold catalyzed matter. The curve corresponding 
to s = 0.5, Yi = 0.4 is unphysical, but has been added in or- 
der to visualize the importance of trapped lepton number at 
subnuclear densities. 
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It may be useful to compare our subnuclear EOS for 

PN^ witl i those used by other authors. The subnuclear 
EOS^o#"RNS, used in the papers of Hashimoto et al. (1995) 
aii^ Bbmbaci et al. (1995, 96) is very different from that 
iisc&Xn the present paper. In particular, Hashimoto et al. 
(1^95) used a cold (T = 0) EOS for the densities below 
g cir On the other hand, Bombaci et al. (1995, 
'96) stop ;heir hot EOS at the edge of the liquid interior, 
and use t he T = (cold catalyzed matter) EOS for the 
densities-below 0.08 fm~^; in this way, they seriously un- 
derestims tc thermal and neutrino trapping effects on the 
radius ofPNS. 

Our E OS above nuclear density are plotted in Fig. lb. 
The prese nce of a trapped lepton number softens the EOS, 
while thermal effects always stiffen it with respect to that 
for cold catalyzed matter. The softening of the supranu- 
clear E03 at fixed Yj is due to the fact, that a signifi- 
cant trap;ped lepton number increases the proton fraction, 
whi' cli imp lies the softening of the nucleon contribution to 
the feoS. 

It should be stressed, that in contrast to Hashimoto et 
al. (1995) and Bombaci et al. (1995, 96) we used a uni- 
fied dense matter model, valid for both supranuclear and 
subnuclear densities. Also, the fact that we use various as- 
sumptions about the T and s profiles within PNS, enables 
us to study the relative importance of the temperature 
profile and that of a trapped lepton number, for the PNS 
models. 

The mass-radius relation for the static PNS models cal- 
culated using various versions of our EOS for the hot inte- 
rior is shown in Fig. 2. We assumed rij^ = 5 x 10~^ fm~^, 
which was consistent with our definition of the "neutri- 
nosphei^' . For the sake of comparison, we show also the 
massr^dius relation for the T = (cold catalyzed matter) 
EQw; wh.ch corresponds to cold neutron star models. In 
twi case of the isothermal hot interior with central tem- 
^^;iperature^rc = 25 MeV we note a very small increase of the 
maximuni mass, as compared to the T = case (c.f., Bom- 
baci et al. 1995, 1996). However, the effect on the mass- 
radius relation is quite strong, and increases rapidly with 
decreasing stellar mass. In the case of the isentropic EOS 
with a triipped lepton number, [s = 2, Yj = 0.4], the soft- 
ening ofthe high-density EOS due to the trapped Yi leads 
to the decrease of M^ax compared to the T = case; as far 
as the valae of M^ax is concerned, the softening effect of Yi 
prevails over that of finite s (this is consistent with results 
of Takatsuka 1995 and Bombaci et al. 1995, 1996). How- 
thermal effect on the radius is very important 
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Fig. 1. b. Pressure versus baryon density for our model 
of dense hot matter, under various physical conditions, for 
supranuclear densities {n > 0.16 fm~^). Notation as in Fig. 
la. 



even in case of = 0.4. This can be seen by compar- 
ing the [s = 2, Yi = 0.4] curve with that corresponding to 
the unphysical, fictitious case of [s = 0.5, Yj = 0.4]. 

The very initial state of a PNS corresponds to a sig- 
nificant trapped lepton number. With our assumption of 
a "standard" composition of dense matter (i,e., exclud- 
ing large amplitude /<'~-condensate, or a very large per- 
centage of hyperons in cold dense matter), the maximum 
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Fig. 2. The gravitational mass versus stellar areal radius for 
static models of the protoneutron stars and neutron stars, 
under various assumptions concerning the physical condi- 
tions within the stellar interior. The curve corresponding to 
s — 0.5, Yi = 0.4 is unphysical, but has been added in order to 
visualize the relative importance of the trapped lepton num- 
ber and thermal effects. The curve T = corresponds to cold 
neutron stars. The curve Tc = 25 MeV includes the factor. 



baryon mass (baryon number) of PNS is lower than that 
of cold NS. Therefore, in our case a stable PNS trans- 
forms into a stable NS, and the scenario PNS — >Black 
Hole, considered by Baumgarte et al. (1996) is excluded. 

At a given mass, the radius of a PNS is significantly 
larger than that of a cold NS. As remarked by Hashimoto 
et al. (1995), this should have important implications for 
rotating PNS. It should be stressed, however, that the 
value of radius, especially for PNS which arc not close to 
the Mmax configuration, turns out to be quite sensitive 
to the location of the edge of the hot neutrino-opaque 
interior (i.e., to the value of n,y). The choice of Bombaci 
et al. (1996) would lead to a much smaller effect on R, 
while that of Hashimoto et al. (1995) would result in larger 
values of the PNS radii. 



2.5. Stationarity 

The EOS of PNS is evolving with time, due mainly to the 
deleptonization process, which changes the composition of 
matter, and also due to changes of the internal temper- 
ature of the star. However, these changes occur on the 
timescales Tevoi ~1-10 s, which are three or more orders of 
magnitude longer than the dynamical timescale, govern- 
ing the readjustment of pressure and gravity forces. This 



dynaniica l timescale Tdyn ~ 1 ms corresponds also to the 
ic periods of the PNS pulsations and of their 
. In view of this, we are able to decouple 
ion from its dynamics, and treat its rotation 
approximation, with a well defined EOS 
S matter. 

the neglected dynamical processes, implied by 

ivc processes and the evolution of the thermal 
of a rotating PNS, is the meridional circulation 
. Strictly uniform rotation is incompatible 
jissumption of a steady thermal state, resulting 
(radiative) equilibrium (see, e.g., Tas- 
). The requirement of radiative equilibrium will 
imply the existence of a meridional circulation 
. However, the velocity of this meridional cir- 
be of the order of the stellar radius divided 
timescale (the timescale of changes of the 
the PNS interior, which is of the order of neu- 
trino difBftion timescale), which is much smaller than the 
rotational velocity (see, eg., Section 8 of Tassoul 1978). In 
view of this, we can neglect the effect of the meridional 
circulation when calculating the mechanical equilibrium 
of a rapidly rotating PNS. Finally, let us notice that in 
a special, idealized case of T* = const, the radiative flux 
vanishes. Then, pure uniform rotation can be realized as 
a steady state of a PNS. 

3. Formulation of the problem 

The problem of the calculation of the stationary state of 
uniform rotation of cold neutron stars, within the frame- 
work of general relativity, was considered by numerous au- 
thors (see the review article of Friedman and Ipser 1992, 
and references therein). Extensive calculations for a broad 
set of realistic equations of state of cold dense matter were 
recently presented in (Cook et al. 1994) and (Salgado et 
al. 1994). Here, wc will extend the methods used at T = 
to the case of hot PNS. We will use the notation and for- 
malism developed in the paper of Bonazzola et al. (1993), 
hereafter referred to as BGSM. 

3.1. Equation of stationary motion 

One of the problems introduced by finite temperature is 
that, if one does not make any other assumption about 
the equilibrium, the equation of stationary motion does 
not have a first integral (Bardeen 1972). 

In the notation of BGSM, the equation of stationary 
motion (Eq. 3.25 of BGSM) reads, when the effects of 
temperature and rigid rotation are included, as 



N 



: Te-"d., 



(4) 



where s is the entropy per baryon, T is the temperature, 
and H = ln[(e + p)/{nmoc'^)] is the so-called pseudo- 
enthalpy (or log-enthalpy), N is the lapse function appear- 
ing in the space-time metric, and T is the Lorentz factor 
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due to rotation. The equation of the stationary motion is 
to be supplemented with the equations determining the 
metric functions (see BSGM for the derivation and the 
explicit form of the complete set of equations). 

In the expression for H, e is the energy density (which 
includes rest energy of matter constituents) and mo is the 
nucleon rest mass. 

This equation is the general relativistic equivalent of a 
well-known Newtonian formula (see e.g. Tassoul 1978). It 
is straightforward to show that a sufficient condition for 
(^ to be integrable is T = T(n). In this case, one obtains 
a first integral of motion of the form : 

N f 

i? + In — - / Te-"ds = const. (5) 

This enables us to calculate the density profile in the en- 
velope of the PNS, because we have assumed a specific 
T{n) profile within it (see section 3.3). Let us stress that, 
due to the integral term in (||), the pseudo-enthalphy is 
no more an explicit function of the metric potentials, as 
it was in BGSM. In fact, one must solve (||) to have H 
within the star. 

In the two particular cases, chosen by us for the hot 
interior, specific first integrals of Eq. (|^) can be found. 
First, in the case of general relativistic thermal equilib- 
rium, T* = Ty = const., it is easy to show that 

N 

fi = II— ~ const., (6) 

where /i = {e+p)/n — Ts is the baryon chemical potential, 
is indeed the integral of (^. The constancy of T* and 
/i* corresponds to the general-relativistic thermodynamic 
equilibrium, if we neglect the time dependence of the EOS 
of PNS (i.e., if we freeze transport phenomena). Second, 
in the case of isentropic profile s{n,T{n)) = const., the 
thermal term in (Q) vanishes and the first integral of Eq. 
(^ reads 

T*s + fJ.* = const. (7) 

3.2. Instabilities and maximum mass 

The instabilities that could develop in rapidly rotating 
PNS, and which could limit the maximum angular fre- 
quency of these objects, are of secular and of dynamic 
type. 

Let us start with the problem of secular stability with 
respect to the axi-symmetric perturbation of rotating con- 
figurations. We denote the total baryon number, total an- 
gular momentum, and total entropy of a PNS by A'^bar, 
J, and S, respectively. We used the secular instability cri- 
terion of Friedman et al. (1988), extended to the case of 
finite temperature. For the purpose of completeness, we 
restate here their lemma (and correct a misprint of their 
paper). 



Lemma : Consider a three-parameter family of uni- 
formly rotating hot stellar models having an equation of 
state of the form p = p{e, T). Suppose that along a conti- 
nous sequence of models labelled by a parameter A, there 
is a point Aq at which iVbar = diVbar/dA, J and S vanish 
and where j^{il*Nh!.r + + f*S) ^ 0. Then the part of 
the sequence for which fi*Nbar + ^J + T*S > is unstable 
for A near Aq. 

Similarly as in the case of Friedman et al. (1988), this 
lemma follows directly from Theorem I of Sorkin (1982), 
with his function S replaced by —M, the i?" quantities 
replaced by iVbar, J, S and the /3" ones by fi*, fi, T*. 
The conditions of the theorem are fulfilled because stellar 
models are configurations for which M is minimized at 
fixed A'bar, J and S, and because difference in M between 
two neighbouring equilibria can be expressed as (Bardeen 
1970) 

dM = Ai*diVbar + f^dJ + T*dS. (8) 

Let us stress here that this last equation requires either 
T* or s to be constant through the whole star, which fixes, 
in each of these two cases, the temperature profile within 
the stellar model. 

We investigated the stability of our models with a spe- 
cific version of this criterion, in which we choose a conti- 
nous sequence of equilibria to be a sequence at fixed iVbar 
and S. The point of the loss of stability is then simply the 
point of extremal J, i.e. : 

where pc is the central density of the star. 

The instability discussed above is a secular one; it will 
develop on the timescale needed to transport the angu- 
lar momentum within the perturbed model, in order to 
decrease the energy of the star while changing its shape 
and structure. However, the timescale of the PNS evolu- 
tion is quite short (seconds), and it is driven by the same 
transport processes involving neutrinos, as those which 
are needed to destabilize the star via the axi-symmetric 
perturbations. In view of this, one might expect that the 
secular instability described above is not efficient in dis- 
rupting the quasi-stationary configuration of rapidly ro- 
tating PNS. 

However, we should remember that the above con- 
siderations apply to the case of infinitesimal perturba- 
tions. Transport processes would be crucial for remov- 
ing the energy barrier separating the initial, secularly un- 
stable configuration, from the dynamically unstable one, 
which would eventually collapse into a black hole. How- 
ever, newly born PNS are expected to be in a highly ex- 
cited state, in which various modes of stellar pulsations 
are excited. One may expect that the energy contained in 
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these pulsations is sufficient to overcome the energy bar- 
rier separating the actual nietastablc, secularly unstable 
state from the dynamically unstable, collapsing one. In 
view of this, we expect that the secularly unstable config- 
urations should be treated like unstable ones. Therefore, 
the critical configuration, given by Eq.(9), will be thus 
considered as the last stable one. 

A rapidly rotating PNS can be also susceptible to other 
types of secular instabilities. The instability with respect 
to the non axisymmetric perturbations can be driven by 
the gravitational radiation reaction (GRR). However, de- 
tailed calculations performed for hot NS suggest, that at 
the temperature exceeding 10^'^ K these instabilities can 
only slightly decrease the maximum rotation frequency of 
uniform rotation, due to the damping effect of the mat- 
ter viscosity (Cutler et al. 1990, Ipser & Lindblom 1991, 
Lindblom 1995, Yoshida & Eriguchi 1995, Zdunik 1996). 
Secular instabilities of rapidly rotating NS with respect 
to the non-axisymmetric "bar" mode were recently in- 
vestigated by Lai & Shapiro (1995) and Bonazzola et al. 
(1995). While Lai & Shapiro (1995) addressed the problem 
of "bar" instability of newly formed NS, they assumed a 
rather unrealistic "ellipsoidal model" , involving only shear 
viscosity as a source of viscous dissipation, and performed 
their calculations within the Newtonian theory of gravity. 
On the other hand, in their rclativistic calculations Bonaz- 
zola et al. (1995) considered only cold NS, and found that 
the secular "bar" instability can set in before the Keple- 
rian (mass shedding) limit is reached only for sufficiently 
stiff EOS of NS matter. The problem of the "bar" instabil- 
ity of PNS, with realistic EOS of the hot interior, will be 
investigated by us in the future. In any case, inclusion of 
possible additional secular instabilities can only decrease 
the maximum rotation frequency of PNS below the val- 
ues obtained using the simplified approach adapted in the 
present paper. 

The fact of the existence of the maximum mass 
of rotating PNS, M^™*^ (and maximum baryon mass, 

-^bar'max) (^^® Scction 5) , puts a well defined and stringent 
limit on which can be reached by PNS. Configurations 
with Mbar > -^bar*max dynamically unstable: no sta- 
tionary solution exists above A'/k''°*1,„„. At the same time, 
the angular frequency of rigid rotation cannot exceed the 
Keplerian value, Qk- These two conditions, combined with 
the instability criterion, expressed in Eq. (9), determine 
the maximum frequency (or, strictly speaking, an upper 
bound on the frequency ) of rigid rotation of PNS. Prac- 
tical implementation of these criteria will be described in 
Section 5. 

3.3. Isothermal and isentropic temperature profiles 

A consistent study of PNS would require an exact treat- 
ment of the thermal transport in the frame of a non- 
stationary spacetime. Unfortunately, such a study is be- 



yond the scope of the present work, and we have thus cho- 
sen to impose "by hand" the temperature profile within 
the star, following prescription described in the subsection 
2.3 . We divided the interior of PNS into the hot interior 
(n > riy), a layer corresponding to the temperature drop 
within the "neutrinosphere" (rij^ — An^, < n < n^), and 
the low temperature, neutrino-transparent outer envelope 
with n < n^ — An^. We have chosen two types of tempera- 
ture profiles within the hot, neutrino-opaque core (n > n^,) 

— the isothermal profile : T = T*jj, with T* = const., 

— the isentropic profile : T = T{n), with s{n,T{n)) = 
const. 

For the transition region and the low temperature en- 
velope, we used a suitable profile T(n) : 

T{n) = 0.2 McV for n < n^ ~ An^ 
T{n) = f{n) for n^ > n > ni, — An^. 

The function / was chosen for computational convenience 
as a suitable combination of an exponential and a gaussian 
function, selected to lead to T{n) of class through the 
transition (temperature drop) region. Our typical choice 
was An^/n^ = 2 10~^; increasing this value up Anyjny = 
0.2 led to a very small increase of the stellar radius. 

As the mass of a massive PNS is almost entirely con- 
tained in the hot neutrino-opaque core (a rough estimate 
for a typical star gives less than 10~^ of the total mass for 
the cool envelope mass), we supposed that our stability 
criterion remained valid also in the case of the presence of 
the low temperature envelope. 

4. Numerical method 

We used a code based on the 3-1-1 formulation of the 
Einstein equations in stationary axisymmetric spacetimes 
(BGSM). The four elliptic equations obtained were solved 
by means of a spectral method, in which the functions are 
expanded in different polynomial bases (Chebyshev for r, 
Legendre for 6 and Fourier for (f). We refer the reader to 
BGSM for a detailed description of the code, including the 
description of the "virial" indicator used for monitoring 
the convergence and the precision reached. Let us just say 
that the code was modified to take into account thermal 
effects and, contrary to Salgado et al. (1994), to converge 
to the solution with various quantities being held fixed (for 
example Mbar or J or 5). Here we briefly outline some of 
the numerical checks we made. 

The two-parameters EOS was interpolated using bicu- 
bic splines subroutine from the NAG library. In this way, 
the thermodynamic functions are of class C^, but the ther- 
modynamic consistency is not conserved. 

The global relative error, evaluated by the means of the 
virial check (see BGSM), is ~ 2-10-3. This relatively "low" 
precision is due to the thermodynamical inconsistencies. 
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We used two grids in r for the star and one in 1/r 

for the exterior with N,. = 33 in each, and one grid in 9 
with N$ = 17. We checked that a greater Nr or Ne does 
not change the results by more than a few 10~^ at most, 
which stays within the global precision reached. Let us 
stress that such a low number of grid points is sufficient 
to reach high accuracy within a spectral method (that 
would not be the case with a finite difference scheme). 

The temperature drop at the "neutrinosphere" is not 
always located on the border of the internal grid in r, 
which could influence the precision (remember that the 
spectral methods are very sensitive to discontinuities). We 
checked that, in fact, even when the "neutrinosphere" is 
far from the border of the grid, the global physical quanti- 
ties of the star did not change much (we found also relative 
variations of a few 10~* at most). 

For a "simple" model where [TcPc^cj were held 
fixed, convergence required no more than 50 iterations 
and ~ 1 minute of CPU time on a Silicon Graphix Indigo^ 
workstation. However, the number of iterations can reach 
400 if one needs to converge to, e.g., a Keplerian configu- 
ration at fixed supra-massive baryon mass. 

Finally, we compared our results at T = and f7 = 
with the results of another code used by one of us (P. 
Haensel), and found that the relative differences in the 
global properties of the stellar models are of the order of ^ 
10-3, which can be imputed to the different interpolation 
procedures. 




12 14 16 n 

Central temperature T^, [MeV] 



Fig. 3. The value of the Keplerian frequency, ^k/'2i^, versus 
central temperature, Tc, for protoneutron star models with an 
isothermal, neutrino opaque, = Q (zero trapped lepton num- 
ber) hot core. The baryon mass of rotating models is fixed at 
A^bar = 1.4 Mq. Different lines correspond to different values 
of the nucleon density at the edge of the hot core, . The value 
of increases by 2 10"^ fm"^ when going upwards from one 
line to the next one. 



5. Numerical results. Maximally rotating pro- 
toneutron stars 

For a given EOS, rapidly rotating PNS models can be 
divided into two families: normal and supramassive one. 
Normal models are those with baryon mass which does 
not exceed the maximum allowable baryon mass of static 
(non-rotating) configurations, -^bar'^max- ^ normal rotat- 
ing model can be transformed into a static configuration 
of the same baryon mass, through a continuous decrease 
of n. Rotation of PNS increases their maximum allowable 
mass with respect to the static case, up to Afbar*max- 

fating PNS models with M^f^^, < M^ar < M^^l^^ are 
called supramassive. Such a supramassive rotating model 
cannot transform into a static model, because in the pro- 
cess of decreasing angular velocity it will collapse into a 
black hole. 

The maximum angular velocity of the normal and 
supramassive rotating models results from different sta- 
bility conditions. In the case of normal rotating config- 
urations the value of f2 is bound by the mass shedding 
limit, which corresponds to the Keplerian velocity at the 
stellar equator, fJ/f (Mbar)- The value of ^Ix turns out to 
be rather sensitive to the location of the "neutrinosphere" 
within the PNS. This sensitivity is particularly large in the 
case of the isothermal profile of the hot neutrino-opaque 



core. This is visualized in Fig. 3, where we show the depen- 
dence of the mass shedding limit il^ for PNS with baryon 
mass 1.4 Mq on the value of for the central temper- 
atures ranging between 10 and 20 MeV. The dependence 
on the value of weakens with increasing Mbar- Also, 
this effect is much less important in the case of isentropic 
PNS. 

For supramassive rotating PNS, the limit on Q. has 

been set by the condition of stability of rotating models 
with respect to the axi-symmetric perturbations, which 
was combined with the mass shedding stability condition. 
Our calculations have been performed for both isothermal 
(T* = const.) and isentropic (s = const.) hot interiors of 
PNS. In both cases, the absolute maximum of Keplerian 
frequency for rotating models, which were stable with re- 
spect to the axi-symmetric perturbations, was obtained for 
a rotating configuration with a maximum baryonic mass 
(and gravitational mass), M^^°^^^^^ [M^l^]. Actually, ro- 
tating configuration with A'/max'' and that with finiax do 
not generally coincide (see, e.g.. Cook et al. 1994, Ster- 
gioulas & Friedman 1995). However, the difference is very 
small, and it could not be detected within the precision 
of our numerical code. The value Mmax"* depends on the 
value of T* in the case of isothermal hot interior, but we 
preferred to parametrize it in terms of central tempera- 
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Table 1. Parameters of the static and rotating maximum mass configurations of protoneutron stars 



EOS 


^ ^-(stat) 
-'^^max 


j-)(stat) 
^max 




p(rot) a 


^max / 27r 




[Me] 


[km] 


[Me] 


[km] 


[Hz] 


T = 0, n = 


2.048 


10.59 


2.430 


14.34 


1625 


Te = 25, = 


2.053 


11.17 


2.322 


14.79 


1521 


s = 0.5, Yi = 0.4 


1.957 


10.85 


2.180 


14.40 


1522 


s = 2,Yi= 0.4 


1.977 


11.53 


2.172 


15.32 


1388 



" Equatorial radius of maximaly rotating configuration. 



ture, Tc = T(r = 0). The dependence of Mmax on the 
vahic of Tc is displayed in Fig. 4 a. The effect of tempera- 
ture on Mmax is opposite to that seen for Mmax*' : thermal 
effects lower the value of Mmax' as compared to that for 
cold neutron stars. This is due to the thermal increase 
of the equatorial radius, which prevails over the thermal 
stiffening of the central core of PNS. The dependence of 
Mmax (rot) on the value of s in the case of isentropic hot 
interior of PNS is represented in Fig. 4 b, in the case of 

7^ (trapped lepton number) and for Yi = 0.4. The 
decrease of Afmaxl^ot) is there smaller than in the case of 
the isothermal PNS models. 

In Fig. 5 a, 5 b we show our results for the maxi- 
mum rotation frequency of stable PNS models, reached 
for supramassive configurations. In the case of the isother- 
mal PNS, thermal effects tend to decrease the value of 
f^max! but even in the case of Tc = 25 MeV the relative 
effect does not exceed ten percent. Irregularities at lower 
Tc result from the limited accuracy of determination of 
the critical configuration with an isothermal core. 

Maximum rotation frequency in the case of the isen- 
tropic hot neutrino-trapped cores of PNS is shown, for 
several values of the central entropy per baryon, in Fig. 5 
b. The value of ^^max is also decreasing with increasing s. 

In general, we find that the decrease of J^max, as com- 
pared to cold neutron stars, is in our case smaller than that 
found by Hashimoto et al. (1995); this may be due to the 
lower value of the density at the edge of the hot core, as- 
sumed by these authors. Some differences may also be due 
to different supranuclear EOS, and to a different treatment 
of the thermal effects, in particular, to their assumption 
of T = const.. In order to visualize the importance of the 
relativistic effects on T{r) in the isothermal interior of a 
PNS, we show in Fig. 6 the temperatures T(r) and T* , 
for Mbar = 1-5 Mq, and a uniform Keplerian rotation. As 
one can see, in the central part of the PNS, T is about 
30 percent greater than T*. A similar effect is seen also 
in the case of an isentropic hot interior, which in Fig. 6 
corresponds to s = 2, Yj = 0.4. 



6. Empirical formula for /2max 

The calculation of the rotating PNS (and NS) models is 

incomparably more difficult than that of the static models. 
In the case of cold NS, one finds a surprisingly precise uni- 
versal formula, which relates the maximum rotation fre- 
quency to the mass and radius of the maximum mass con- 
figuration for the static models (Haensel & Zdunik 1989, 
Friedman at al. 1989, Shapiro et al. 1989, Haensel et al. 
1995, Nozawa et al. 1996), 

/ ,^(stat) \ 5 / p(stat) \ ^1 

""" = ^(^] [W^) ■ W 

where the most recent value of C, based on calculations 
performed for a broad set of realistic cold EOS of dense 
matter, is Ccoid = 7750 s"^ (Haensel et al. 1995). The 
"empirical formula", Eq. (10), reproduces the values of 
f^max for cold NS models with typical precision better than 
5% (although in some specific cases the deviations can 
reach nearly 7%, see Nozawa et al. 1996). 

The validity of the empirical formula for fimax is of 
great practical importance. For cold EOS of dense matter, 
it enables one to get immediately a rather precise estimate 
of the value of f2max, using easily calculated static neutron 
star models, and avoiding in this way incomparably more 
difficult 2-D calculations of rotating NS models and the 
analysis of their stability. 

It is of interest to check, whether "empirical formula" 
is valid also in the case of hot PNS. Let us notice, that 
the subnuclear EOS of PNS, with possible lepton num- 
ber trapping, is very different from that of cold catalyzed 
matter (see Fig. 1). In order to check the validity of the 
empirical formula, we define the EOS depending parame- 
ter Ceos, 

/ , -^(stat) \ ~i / n(stat) \ f 

where the right-hand-side is calculated using exact results 
for a specific EOS. The values of Ceos for our models of 
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Fig. 4. a. Maximum gravitational mass of rotating pro- 
toneutron star models with isothermal cores and = 
(zero trapped lepton number), versus central temperature, 
Tc. The nucloon density at the outer edge of the hot core 

= 5 10"-^ fm^^ 



Fig. 5. a. Maximum rotation frequency of stable protoneu- 
tron star models with isothermal hot cores (with zero trapped 
lepton number), versus central temperature, Tc. The nucleon 
density at the outer edge of the hot core is rii/ = 5 10~^ fm~^. 
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Fig. 4. b. Maximum gravitational mass of rotating protoneu- 
tron star models with iscntropic hot cores with Yi — 0.4, versus 
central entropy, s. The nucleon density at the outer edge of the 
hot core is ni/ = 5 10~^ fm~*. 



Fig. 5. b. Maximum rotation frequency of stable protoneutron 
star models with isentropic hot cores with Y; = 0.4, versus 
central entropy, s. The nucleon density at the outer edge of 
the hot core is rii, = 5 10~^ &n~^. 
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Fig. 6. Difference between T and T* in the hot interior of a 
A<fbar = 1-5 Mq PNS iu Kcplcrian rotation, for the two models 
of the hot PNS core. Dashed lines correspond to the isothermal 
(Tc = 15 MeV), and the dash-dot lines to the isentropic model, 
respectively. For each pair of curves, upper one corresponds to 
T(r), and the lower one to T*{r). Temperatures are plotted 
versus radial coordinate r in the equatorial plane of rotating 
PNS. 



PNS are displayed in Fig. 7 a, 7 b. Let us notice, that 
for isentropic hot cores with trapped lepton number the 
vahie of Ceos shows a decreasing trend with increasing 
s. In the case of isothermal hot cores with Y^jy = we see 
the opposite trend: the value of Ceos tends to increase 
with central temperature. The value of Choi = 7800 s^^ 
will lead to a very good empirical formula for PNS, with 
precison comparable to that used for cold NS. Within the 
precision of the empirical formulae, the values of Ccoid and 
Chot can be considered as being equal. 



7. Numerical results. Prom hot protoneutron stars 
to cold neutron stars 



Rapidly rotating hot, neutrino-opaque PNS evolves even- 
tually into a cold, rotating NS, which under favourable 
circumstances can be observed as a solitary pulsar. Let us 
assume, that such a transformation took place at constant 
baryon mass: this assumption is valid if mass accreted af- 
ter the formation of a PNS is negligibly small. (Notice, 
that we can define the moment of formation of a PNS 
as that at which the accretion ends). Second assumption 
refers to the angular momentum of the star. Let us ne- 
glect for the time being the angular momentum loss due 




Central temperature 



Fig. 7. a. The cocfRciont Ceos, given by Eq.(ll), versus cen- 
tral temperature, Tc, for isothermal models of the zero trapped 
lepton number (y„ = 0) hot cores of protoneutron stars. 




1 1.5 2 

Central entropy per baryon s [kg] 



Fig. 7. b. The coefficient Ceos, given by Eq.(ll), versus cen- 
tral entropy per nucleon, for isentropic Yj = 0.4 models of hot 
cores of protoneutron stars. 
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Table 2. Maximum rotation frequencies and maximum baryon 
masses of PNS with various EOS, and corresponding frequen- 
cies of cold NS configurations, obtained via cooling at constant 
J and Mbar (J-losses due to neutrino emission are neglected). 



EOS 


bar, max 

[Mq] 


[Hz] 


^max /27r 

[Hz] 


[Hz] 










t; 


Te = 9, = 


2.806 


1603 


1582 


\ 


Te = 25, n = 


2.695 


1521 


1403 


a 


s = 2, y^ = 0.4 


2.172 


1388 


1112 





to emission of neutrinos (mainly during deleptonization) , 

as well as that due to gravitational radiation (if rotating 
star deviated from the axial symmetry). Inclusion of both 
effects could only decrease the final angular momentum of 
the NS. Under these two assumptions, the transformation 
of a PNS into a solitary radio pulsar takes place at con- 
stant baryon mass, Mbar, and angular momentum, J. The 
maximiim angular frequency of a solitary radio pulsar of 
a given baryon mass would be then determined by the pa- 
rameters of a maximally rotating PNS of the same baryon 
mass. 

In the present paper we restrict ourselves, similarly 
as Hashimoto et al. (1995), to the case of rigid rotation 
of PNS. We will first consider "typical" range of baryon 
masses expected for the NS born in a gravitational col- 
lapse of a massive stellar core, 1.4 Mq < Mbar < 2.0 Mq 
(for our EOS of dense matter this would correspond to 
1.3 M0<M< 1.8 Mq). All stars considered are normal, 
and the maximum initial frequency of PNS will be then 
f7A'[hot,y^ = 0], and f}if[bot,yi = 0.4], for the limit- 
ing case of the zero trapped lepton number, and maxi- 
mum trapped lepton number neutrino-opaque cores, re- 
spectively. The values of these limiting angular frequen- 
cies, for the isothermal and isentropic PNS cores, are plot- 
ted versus Mbar in Fig. 8 (thin dashed and dash-dotted 
lines, respectively). 

Our results displayed in Fig. 8 lead to some conclu- 
sions, relevant to those neutron stars which remained soli- 
tary after their birth in the gravitational collapse of a 
massive stellar core. If the starting configuration (at the 
end of significant accretion, after the revival of a stag- 
nated shock) was that with trapped lepton number and 
isentropic, then the maximum frequency of such a neu- 
tron star, in the gravitational mass range of 1.3 - 1.8 Mq, 
cannot exceed 600 - 900 Hz, the lower freqiiency limit cor- 
responding to the lower mass limit. The corresponding 
range of minimum periods would be from 1.6 ms to 1.1 ms, 
respectively. 

If the "initial isolated configuration" (i.e., that with 
no subsequent accretion) was a little older, and therefore 
deleptonized, and if we assume that it was isothermal. 
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Fig. 8. Maximum rotation frequencies versus baryon mass 
of rotating configurations. Solid line: maximum rotation fre- 
quencies of cold neutron star models. Ox (cold). Thin dsished 
line: maximum rotation frequency for isothermal models with 
Tc — 33 MeV with no trapped lepton number {Y^, = 0), de- 
noted in the text as Slicfhot, Yv ~ 0]. Heavy dashed line: angu- 
lar frequency of cold configurations, resulting from those cor- 
responding to thin dashed line by cooling at constant J and 
Mbar. Thin dash-dotted line: maximum rotation frequency for 
isentropic models, denoted in the text as Q,k[s — 2,Yi = 0.4]. 
Heavy dash-dotted line: angular frequency of cold configura- 
tions, resulting from those corresponding to thin dash-dotted 
line by cooling at constant J and Mbar (J-losses due to neu- 
trino emission are neglected). 



then the resulting constraints would be less stringent: for 
the range of masses 1.3 - 1.8 Mq we get minimum rotation 
periods of 1.25 - 1.1 ms, respectively. 

In both cases cooling resulted in speeding up of the ro- 
tation of the star: such situation is characteristic of "nor- 
mal rotating models" . 

Similar evolutionary considerations can be applied to 
the "maximally rotating configurations" : isentropic one 
with s = 2, Yi = 0.4 and isothermal ones with = 
9 MeV, 25 MeV, respectively. Our results are shown in Ta- 
ble 2. The initial "hot" configurations are supramassive, 
and evolve into more compact cool ones (at constant Mbar 
and J) decreasing their angular frequency (slowing down). 
This purely relativistic effect was recently pointed out by 
Hashimoto et al. (1995). The "relativistic slowing down" 
during cooling of the maximally rotating isentropic config- 
uration with s = 2, Yi = 0.4 corresponds to the increase 
of the period from Fmin[s = 2,y; = 0.4] = 0.72 ms up to 
^min = 0-90 "^s- If such a scenario of formation of solitary 
pulsars is valid, an absolute limit on their period would 



14 



J.O. Goussard et al.: Rapid uniform rotation of protoneutron stars 



be -P^irn not Pinin{T = 0) = 0.61 ms, obtained for the 
T = EOS. The slowing down factor of ~ | coincides with 
that obtained by Hashimoto et al. (1995) for a different 
EOS and a different scenario of formation of solitary cold 
pulsars. Actually, they assumed that the initial configura- 
tion is that with T = const. ,Yi, = 0. Their initial model 
would to some extent correspond to our isothermal models 
(remember however the factor ^ o^^ir temperature pro- 
files, which is absent in their calculations). Using Tables 
1, 2 we obtain PUATc = 25, = 0]/P„,i„[T = 0] = 0.92, 
which is quite close to 1. 

8. Discussion and conclusions 

In the present paper we studied rapid rotation of pro- 
toneutron stars, using a specific model of hot, dense mat- 
ter. Our models of rapidly rotating protoneutron stars 
were based on several simplifications, which were neces- 
sary in order to make the problem tractable. In order 
to avoid difficulties and/or ambiguities of the largely un- 
known "real situation" , we restricted ourselves to study- 
ing idealized, limiting cases. In many places we introduced 
approximations, which were crucial for making numerical 
calculations feasible. 

We assumed rigid rotation within protoneutron star. 
This simplified greatly our considerations, and reduced 
dramatically the number of stellar models. Actually, our 
formalism allows calculation of the quasi-stationary differ- 
entially rotating configurations, after introducing an addi- 
tional FdiVl term appearing on the right hand side of the 
equation of stationary motion, Eq. (4) (see Section 5.1 of 
BGSM). We plan to perform studies of differentially rotat- 
ing protoneutron stars in the near future, as the next step 
in our investigations of dynamics of protoneutron stars. 
While existing numerical simulations of gravitational col- 
lapse of rotating cores of massive stars yield a differentially 
rotating protoneutron star, the calcidations stop too early 
after bounce. In contrast to our models of protoneutron 
stars, the object which comes out from the simulations 
of Mocnchmeyer and Mueller has an extended, very hot 
envelope, produced by the shock wave immediately after 
bounce (Janka & Mocnchmeyer 1989a, b, Mocnchmeyer & 
Mueller 1989). Also, the lack of knowledge of the initial 
angular velocity distribution within the collapsing core re- 
sults in the uncertainty in the rotational state of produced 
protoneutron star. Clearly, the study of quasistationary 
differential rotation of protoneutron stars should take into 
account all these uncertainties. 

Deleptonization of protoneutron star is connected with 
energy and angular momentum losses. However, angular 
momentum taken away by neutrinos is expected to con- 
stitute at most a few percent of the total stellar angular 
momentum (Kazanas 1977). Inclusion of this effect would 
slightly decrease some of our "evolutionary limits" of Sec- 
tion 7. 



Our treatment of the thermal state of the protoneu- 
tron star interior should be considered as very crude. The 
temperature profile might be affected by convection. Also, 
our method of locating the "neutrinosphere" was very ap- 
proximate. Clearly, the treatment of thermal effects can 
be refined, but we do not think this will change our main 
results. 

Our calculations were performed for only one model of 
the nucleon component of dense hot matter. The model 
was realistic, and enabled us to treat in a unified way the 
whole interior (core as well as the envelope) of the pro- 
toneutron star. However, in view of the uncertainties in the 
EOS of dense matter at supranuclear densities, one should 
of course study the whole range of theoretical possibilities, 
for a broad set - from soft to stiff - of supranuclear, high 
temperature EOS. An example of such an investigation, 
in the case of the static protoneutron stars, is the study of 
Bombaci et al. (1996). In view of the possible importance 
of the protoneutron star - neutron star connection for the 
properties of solitary pulsars, similar studies should also 
be done for rotating proto-neutron stars. 

The calculations of the present paper were done under 
the assumption, that the frequency of uniform rotation 
of protoneutron stars is limited only by the mass shed- 
ding and the secular instability with respect to the axi- 
symmetric perturbations. In view of this, our results can 
be considered only as upper bounds to the maximum rota- 
tion frequency. Our main conclusion is that the minimum 
rotation period of solitary neutron stars, born as rapidly 
rotating protoneutron stars, is significantly larger, than 
the corresponding limit for cold neutron stars. Inclusion 
of additional secular instabilities in rapidly, uniformly ro- 
tating protoneutron stars can only strengthen this conclu- 
sion. 
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